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Abstract 

We propose to use the hadron number fluctuations in the hmited momentum regions to study the 
evolution of initial flows in high energy nuclear collisions. In this method by a proper preparation 
of a collision sample the projectile and target initial flows are marked in fluctuations in the number 
of colliding nucleons. We discuss three limiting cases of the evolution of flows, transparency, 
mixing and reflection, and present for them quantitative predictions obtained within several models. 
Finally, we apply the method to the NA49 results on fluctuations of the negatively charged hadron 
multiplicity in Pb+Pb interactions at 158A GeV and conclude that the data favor a hydrodynamical 
model with a significant degree of mixing of the initial flows at the early stage of collisions. 
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1. The main goal of investigations of high energy nucleus-nucleus (A+A) collisions is to 
uncover properties of strongly interacting matter at high energy densities and, in particular, 
to look for its hypothetical phases and transitions between them. Qualitative features of 
the rich experimental data collected thus far indicate that the produced matter experiences 
strong collective expansion and it is close to local equilibrium Moreover, the properties 
of the matter change rapidly at the low CERN SPS energies {^/sNN ~ 8 GeV) suggesting 
the onset of deconfinement and thus the existence of a new state of matter, a Quark Gluon 
Plasma . The properties of this new phase are under active studies in A+A collisions 
at the BNL RHIC ^ = 200 GeV). 

We are, however, far from a full understanding of the A+A dynamics. Many models 
based on different assumptions compete with each other and a consistent description of all 
aspects of the data within a single model is missing. The largest uncertainties concern the 
early stage of collisions. It is unclear how initial nuclear flows of energy and charges evolve. 
The majority of the dynamical models (e.g. the string-hadron transport approaches and the 
quark-gluon cascade models 0,0,01) predict or assume that the colliding nuclear matter is 
transparent. The final longitudinal flows of the hadron production sources or the net baryon 
number related to the projectile and target follow the directions of the projectile and target, 
respectively. We call this class of models transparency (T-)models. Since the pioneering 
works of Fermi and Landau statistical and hydrodynamical approaches are successfully 
used to describe high energy nuclear collisions. Many models within this group, including the 
first Fermi formulation, assume full equilibration of the matter at the early stage of collisions. 
The initial projectile and target flows of energy and charges are mixed. The approaches which 
predict or suppose the full mixing of the projectile and target flows we call the mixing (M- 
) models. Let us note that there are models which assume the mixing of hadron production 
sources (inelastic energyVwhereas the transparency of baryon number flows, e.g. statistical 
model of the early stage ^ and the three- fluid hydrodynamical model Q]. Finally, one may 
even speculate that the initial flows are reflected in the collision process, i.e. the flows of 
matter related to the target and the projectile change their directions. This class of models 
we call the reflection (R-)models. The sketch of the rapidity distributions resulting from the 
T-, M- and R-models are shown in Fig.^ The spectra related to the projectile and the target 
can be easily distinguished in the figure because they are marked in color and hatching the 
same way as the initial projectile and target nuclei. In this paper we propose a method to 
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FIG. 1: Sketch of the rapidity distributions of the baryon number or the particle production sources 
(horizontal rectangles) in nucleus-nucleus collisions resulting from the transparency, mixing and 
reflection models. The spectator nucleons are indicated by the vertical rectangles. In the collisions 
with the fixed number of projectile spectators only matter related to the target shows significant 
fluctuations (vertical arrows). 

mark the matter related to the projectile and the target in fluctuations (the MinF- method), 
which allows to test experimentally different scenarios of the collision process. Finally we 
apply the MinF-method to the NA49 experimental data on Pb+Pb collisions at 158A GeV 
(v^ = 17.2 GeV) [nl. 

2. In each A+A collision only a part of all 2 A nucleons interact. These are called partic- 
ipant nucleons and they are denoted as N^"-^ and N^p^^ for the projectile and target nuclei, 
respectively. The nucleons which do not interact are called the projectile and target spec- 
tators, Ng^"-' = A — N^"-' and A^^"''^ = A — N^p"^^ . The fluctuations in high energy A+A 
collisions are dominated by a trivial geometrical variation of the impact parameter. How- 
ever, even for the fixed impact parameter the number of participants, Np = N^"-' + N^p"^ , 
fluctuates from event to event. This is caused by the fluctuations of the initial states of the 
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colliding nuclei and the probabilistic character of an interaction process. The fluctuations 
of Np usually form a large and uninteresting background. In order to minimize its contribu- 
tion NA49 selected samples of collisions with flxed numbers of projectile participants. This 
selection is possible due to the measurement of Ng^"-^ in each individual collision by use of a 
calorimeter which covers the projectile fragmentation domain. However, even in the samples 
with Np'°-' = const the number of target participants fluctuates considerably. Hence, an 
asymmetry between projectile and target participants is introduced, i.e. N^°-' is constant, 
whereas N^p^^ fluctuates. This difference is used in the MinF-method to distinguish between 
the flnal state flows related to the projectile and the target. Qualitatively, one expects large 
fluctuations of any extensive quantity (e.g. net baryon number and multiplicity of hadron 
production sources) in the domain related to the target and small fluctuations in the projec- 
tile region. When both flows are mixed intermediate fluctuations are predicted. The whole 
procedure is presented in a graphical form in Fig. ^ Clearly, the fluctuations measured in 
the target momentum hemisphere are larger than those measured in the projectile hemi- 
sphere in T-models. The opposite relation is predicted for R-models, whereas for M-models 
the fluctuations in the projectile and target hemispheres are the same. 

This general qualitative idea is further on illustrated by quantitative calculations per- 
formed within several models ordered by an increasing complexity. 

3. Let us begin with considering fluctuations of the net baryon number measured in differ- 
ent regions of the participant domain in collisions of two identical nuclei. These fluctuations 
are most closely related to the fluctuations of the number of participant nucleons because 
of the baryon number conservation. In the following the variance, Var{x) = (x^) — (x)^, 
and the scaled variance, = Var{x)/ (x), where x stands for a given random variable 
and (■■■) for event-by-event averaging, will be used to quantify fluctuations. We denote 
by cjp'"^ = Var(N^p^^)/ {N^p'^^) the scaled variance of the number of target participants and 
by ujb = Var{B) / (B) the scaled variance of the net baryon number, B. In each event we 
subtract the nucleon spectators when counting the number of baryons. The net baryon 
number, B = Np — N^, equals then the number of participants Np = N^p^^ + N^"-' . At 
fixed N^"'' , the Np number fluctuates due to the fluctuations of N^p^^ . The distribution in 
N^p^^ can be characterized by its mean value, {N*p^'^) ~ N^"\ and a scaled variance, ufp^^ . 
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Thus, for the net baryon baryon number B one finds, 

Var{Np) {{N'rY) - {N^r 1 targ 

- " — mr^) — - 2 ' 

for the fiuctuations in the full phase space of participant nucleons. A factor 1/2 in the right 
hand side of Eq. ((T)) appears because only a half of the total number of participants fiuctuates. 
Let us introduce uj^ and cj^, where the superscripts p and t mark quantities measured in the 
projectile and target momentum hemispheres, respectively. By assumption, the mixing of 
the projectile and target participants is absent in T- and R-models. Therefore, in T-models, 
the net baryon number in the projectile hemisphere equals N^^"^ and does not fiuctuate, 
i.e. a;^(T) = 0, whereas the net baryon number in the target hemisphere equals Np"-"^^ and 
fiuctuates with ci;^(T) = uj^p^^. These relations are reversed in R-models. 

We introduce now a random mixing of baryons between the projectile and target hemi- 
spheres. Let a be a probability for (projectile) target participant to be detected in the 
(target) projectile hemisphere. It is easy to show that: 

uj'b = {1 - a)2 uo'p + 2a{l - a) , oo^ = + 2a(l - a) . (2) 

A (complete) mixing of the projectile and target participants is assumed in M-models. Thus 
each participant nucleon with equal probability, a = 1/2, can be found either in the target 
or in the projectile hemispheres. In M-models the fiuctuations in both hemispheres are 
identical. The limiting cases, a = and a = 1 of Eq. (j21) correspond to T- and R-models, 
respectively. In summary the scaled variances of the net baryon number fiuctuations in the 
projectile, u^, and target, tu^, hemispheres are: 

uUT) = , uUT) = , (3) 

a;^(M; rr) = c.^(M; rr) = \ + \ , (4) 

o;^(i?) = , u'^{R) = , (5) 

in the T- Q, M- ^ and R- (0) models of the baryon number fiow. When deriving Eq. Q 
we assumed that the baryons are distributed randomly in the rapidity space thus the abbre- 
viation rr in the left hand side of Eq. (@]) stands for random rapidities. This implies that 
even for a fixed number of N^p"^^ , i.e. for u/p^^ = 0, the baryon number in the projectile and 
target hemispheres fiuctuates, a;^(M; rr) = a;^(M; rr) = 1/2. 



In a mixing model in which baryon rapidities do not fluctuate from collision to collision, 
but their positions are fixed (the fixed rapidity, {fr), model) the scaled variances in the 
projectile and target hemispheres read: 

uj^M; fr) = co^M; fr) = ^ . (6) 

Note that the term 1 /2 of the right hand side of Eq. (jU is absent in ® . 

Eventually, for an estimate of the magnitude of the expected fluctuations in different type 
of models we consider an example of Pb+Pb collisions at 158A GeV. The scaled variance 
of the number of target participants at the fixed number of projectile participants (i.e. 
0^°-' = 0) can be calculated within the string-hadronic models. The corresponding results 



12| obtained for the HSD 5] model are shown in Fig. |2l left. Using Eqs. (jHE)) and the 
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FIG. 2: Left: The scaled variance ujp^^ for the fluctuations of target participants Np^^ as a 
function of Np°-' calculated within the HSD model. Right: The scaled variances lo^ versus 
j^proj Qi^i^g^jj^gfi within T- (dashed line) and R-models (dashed-dotted line), Eqs. (jSHSJ. The upper 
solid line shows predictions of the M-model with random rapidities of baryons Q whereas the lower 
solid line corresponds to the M- models with fixed rapidities of baryons © . For tOp the predictions 
for T- and R-models should be exchanged and the lines for the M- models remain unchanged. 

dependence of uu'^p^^ on N^°-' calculated within the HSD model (Fig. |21 (left)), quantitative 
predictions concerning the baryon number fluctuations for different models can be obtained. 
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The resulting dependencies of the scaled variance of the baryon number in the projectile 
hemisphere on N'^°-' are shown in Fig. |21 (right). As expected large fluctuations are seen 
in R-models, intermediate in M-models and there are no fluctuations in T-models. In the 
M-model the scaled variance increases by 1/2 when baryon positions in rapidity are assumed 
to fluctuate. 

4. The T-, M- and R-models for the baryonic flows give indeed very different predictions 
for cu^ and cu^ for the samples of events with fixed values of Np"-'. However, they may be 
difficult to test experimentally as an identification of protons and a measurement of neutrons 
in a large acceptance in a single event is difficult. Measurements of charged particle multi- 
plicity in a large acceptance can be performed using the existing detectors. In particular, 
the first results on multiplicity fluctuations of negatively charged hadrons, A^_, as a function 
of N^p°^ were recently obtained by NA49 Q for Pb+Pb coUisions at 158A GeV. Note that 
at the CERN SPS and lower energies negatively charged hadrons are predominantly (more 
than 90%) tt" mesons. In the following we consider T-, M- and R- scenarios within sev- 
eral approaches to particle production in high energy nuclear collisions. We suppose that a 
part of the initial projectile and target energy, the inelastic energy, is converted into hadron 
sources. Further on, the numbers of projectile and target related sources are taken to be 
proportional to the number of projectile and target participant nucleons, respectively. The 



physical meaning of a partic 
are wounded nucleons (see 



e source depends upon the model under consideration, examples 
3, 13), strings and resonances (see 0,0]), and volume cells of 



the expanding matter at the freeze-out in the hydrodynamical models. For the independent 



sources one finds regarding the scaled variance of i-th particle species (see e.g. 16|): 



oji = u* + {n*) , (7) 

where uj* denotes the scaled variance for i-th hadron species (e.g., i may correspond to h^) 
from a single source, (n*) is the average multiplicity from a single source, and Q^, is the 
scaled variance for the fluctuation of the number of sources. Assuming that the number of 
hadron sources is proportional to the number of participating nucleons, N^, = const ■ Np, 
one gets: 

{^i ) fi* = fi* = UJp = TliUJp , (8) 
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where Ui is the average multiphcity of the i-th species per participating nucleon. Thus the 
scaled variance ((Tj) of the particle number multiplicity in the full phase space is: 

Ui = u* + Up''^ . (9) 

Consequently, the scaled variances of the z-th hadron multiplicity distribution in T-, M- and 
R-models read: 

u^iT) = u* , uliT) = u* + c^r^ , (10) 
uf{M; rr) = u\{M- rr) = u* + + \ uf' ) , (H) 

u\{M- fr) = u\{M- fr) = u* + uf' , (12) 

ufiR) = u* + u';^^ , uliR) = u: . (13) 

Again two different versions of mixing with random rapidities (jllj) and fixed rapidities fll2|l 
of the source positions are possible. As an example of M-models with fixed rapidities of 
the sources let us consider a model which assumes a global equilibration of the matter at 
the early stage of collisions followed by a hydrodynamical expansion and freeze-out. In this 
case particle production sources can be identified with the volume cells of the expanding 
matter at the freeze-out. They can be treated as uncorrelated provided the effects of global 
energy-momentum conservation laws can be neglected. Due to assumed global equilibration 
of the projectile and target flows the fluctuations in the projectile and target hemispheres 
are identical. The model belongs to the class of M-models. In this model there is one to 
one correspondence between space-time positions and rapidities of the hydrodynamic cells. 
Thus, the source rapidities do not fluctuate and the scaled variances of hadrons in the 
projectile and target hemispheres have the form ()12|1 . 

Note that Eqs. (|Tn|l and (|T!^ are strictly valid provided that a source produces particles 
only in its hemisphere. Due of the finite width of the rapidity distribution resulting from the 
decay of a single source this condition is expected to be violated at least close to midrapidity. 
Thus, in order to be able to neglect the cross-talk of particles between the projectile and 
target hemispheres the width of the rapidity distribution of a single source, Ay^,, should be 
much smaller than the total width of the rapidity distribution. 

Some comments concerning Eqs. are appropriate. There is a general similarity of 

the expressions for produced particles and the corresponding expressions for baryons, Eqs. (0- 
E}. There are, however, two important differences. A single source produces particles in a 
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probabilistic way with an average multiplicity (n*) and a scaled variance u*. Consequently, 



it leads to an additional term, a;*, in all expressions for cjf'*, and an additional factor, n. 



appears in terms related to the fluctuations of the number of sources. Following Eq. (jH)) the 
source number fluctuations can be substituted by u/p^^ , and an average multiplicity, (n*), of a 
single source can be then transformed into an average multiplicity per participating nucleon, 
rii. The term, 1/2, in the r.h.s. of Eq. (fTTj) . as that in Eq. is due to the random rapidity 
positions of the sources in M-models. In the hydrodynamical model particle production 
sources can be identified with the volume cells of the expanding matter at the freeze-out. 
The source rapidities do not fluctuate and Eq. (fTT|) is transformed then into Eq. (fT^ . 

We turn now to a discussion of multiplicity fluctuations of negatively charged hadrons 
in Pb+Pb collisions at 158/1 GeV. The value of {NJ)/{Np) = n_ ~ 2 was measured for 
the studied reactions For simplicity we assume u;* ~ 1, this is valid for the poissonian 
negatively charged particle multiplicity distribution from a single source. Note, in p+p 
interactions at SPS energies and in the limited acceptance of NA49 the measured distribution 
is in fact close to the Poisson one ^]|. This gives: 

ujI(T) = cui(i?) ^ 1 , (14) 

a;P(M; rr) = uj'_{M; rr) ~ 2 + ^ cu^"'"^ , (15) 

ujP_{M; fr) = uj'_{M; /r) ~ 1 + ^ u'P (16) 

uj'_{T) = cu^(i?) ~ 1 + 2 . (17) 

The dependence of on N^"-' from Eqs. ()14llT7jl for T-, M- and R-models is presented in 
Fig. El • 

5. In a recent analysis [3| of d+Au interactions at ^/SNN = 200 GeV ^1^ within the 
wounded nucleon model (WNM) [l^ it was found that the wounded nucleon sources emit 
particles in a very broad rapidity interval which results in the mixing of particles from the 
target and projectile sources. In the following we consider predictions of this model with 
respect to multiplicity fluctuations. 

The WNM assumes that i-th particle rapidity distribution in A+A collisions is presented 

as 

^ = iVr^ F^iy) + N'r Ffiy) , (18) 
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FIG. 3: The dependence of the scaled variance of negatively charged particle multiplicity in Pb+Pb 
collisions at 158A GeV on the number of projectile participants, N^°\ in the projectile hemisphere. 
The predictions for T- H14() and R- ((T7I) models are shown by dashed and dashed-dotted lines, 
respectively. The upper solid line, mixing(rr), corresponds to the M-models with random rapidity 
positions of the sources (|15|) . The lower solid line, denoted as mixing(fr), corresponds to the M- 
models with fixed rapidity positions of the sources (jlbf) . For the target hemisphere the lines of T- 
and R-models should be interchanged, whereas the lines of M-models remain unchanged. We take 
u/p'^^ from Fig. (left) for all types of models. 

where Fliy) and Ffijj) are the contributions from a single wounded nucleon (identified with 
a particle source) of the target and projectile, respectively. The model also requires 

F!{y) = Fl{~y) (19) 

in the center of mass system of the collision. The mean number of particles in the rapidity 
interval Ay for collisions with N'p"-' and A'^^"'^^ is given by 

N,{Ay) = N'r I dy F^{y) + I dy Ff{y) . (20) 

J Ay J Ay 



reflection 



mixing(rr) 
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For interaction of identical heavy ions {N'^'^^) ~ N^°-' , and then Eq. ()20|) yields: 

m/^y)) = Nr I dy [Fl{y) + Ff{y)\ . (21) 

Let us consider now fluctuations of Ni{Ay) at a fixed N^°-' . A contribution to the scaled 
variance of Ni{Ay) (plj) due to the fluctuations of N^p'^^ reads: 

1 2 



LydyFliy) 



Var[Ni{Ay)] ^ 

{N,{Ay)) Ldy [F!{y) + F[{y)] 



u'P ^ n\{Ay) a\{Ay) , (22) 



r f dv F^iy) 



/Ay 

where 

lAyd' 

As previously, for simplicity we assume that a single source emits particles according to the 
Poisson distribution, u*{Ay) = 1. This leads to a general expression on the scaled variance 
for a particle of i-th type: 

u.iAy) = 1 + nliAy) aliAy) co';^^ . (24) 

The parameter a* quantifies the amount of mixing of the projectile and target contributions 
and can vary between and 1. For full acceptance. Ay = [—Ymax, Ymax], Eq- transforms 
to Eq. ©. 

The T-, M- and R- limits of the WNM can be formulated in terms of the distribution 
functions of the single nucleon, 

Ff{y-T) = T,{y) e{y) , F^{y;T) = T,{-y) e{-y) , (25) 
F[{y;M) = F^{y;M) = M,(y) , (26) 
F[{y;R) = R,{-y) 6{-y) , Fj = R,{y) 6{y) . (27) 

The scaled variances in the projectile and target hemispheres, cjf = uJi{y > 0) and 
ujj = uji{y < 0), can be found using Eq. (j2H). It follows that in T- and R- models, the tuf'* 
coincide with those given by Eq. (jiup and Eq. (jl3j) . respectively. For the M- models this 
gives the following result: 

cjf (M; WNM) = uo\{M; WNM) = 1 + ^ Jp^ , (28) 
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which is identical to Eq. (fT^ . In the WNM all projectile (target) sources are assumed to 
be identical and their positions are the same and fixed. Therefore, similar to the hydro- 
dynamical model, the term, 1/2, in the r.h.s. of Eq. (jllj) is absent in Eq. (j28p . This is 
the M-model with the fixed rapidity positions of the sources. The mixing in the considered 
version of WNM results from a broad distribution of particles produced by a single source. 
A complete mixing in the WNM means according to Eq. ()26|) that projectile and target 
source functions become identical. 

6. Let us consider fiuctuations in limited phase-space domains in which only fractions, 
g^'*, of all particles in the projectile or target hemispheres are accepted. Then the scaled 
variances in the acceptance, u^^^ and cu*^^, will be different from the and tu*. We start 
with the scaled variances and tu^ of the net baryon number fluctuations. Assuming that 
inside the proiectile and target hemispheres the baryon rapidities are not correlated one gets 
(see e.g.,|l6|,y): 

= 1 - + ■ ^^'^ • (29) 

It can be shown that the scaled variance of the produced particles in the limited momen- 
tum acceptance within the M-model with flxed source rapidities f|T^ reads: 

u^^,aUM, ^) = l + ^^^4"'■^ (30) 

where {Ni)acc is a mean multiplicity of a particle of i-th type in the acceptance. The formula 
above assumes that the produced particles are uncorrelated in the momentum space, i.e. it 
neglects effects of motional conservation laws and resonance decays. The scaled variance in 
a limited rapidity acceptance. Ay, within WNM can be directly obtained from Eq. ()24j) . It 
coincides with that of Eq. (j30|l . 

Let us consider now as an example the NA49 acceptance, which is located in the projectile 
hemisphere about one and half rapidity units from mid-rapidity, Ay = [1.1; 2.6] in the c.m.s. 
The acceptance probability was measured to be q'^ ^ 0.4 |ll|,ll7[ (i.e. about 40% of negatively 
charged particles in the projectile hemisphere are accepted). In the limiting case of the fixed 
rapidities of the sources, this is assumed to be valid for both the hydrodynamical and WN 
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models, one finds: 



<acc{T) ^ 1, (31) 

c.^,_(M; /r) ^ 1 + 0.2 u;^^ , (32) 
{R; /r) ~ 1 + 0.8 uo'P . (33) 



p 



The relations, n^_{Ay;T) = 0; n^(A?/;M) = 0.5 n_; n^_{Ay; R) = qP n_, with 
~ 0.4 and ?T_ ~ 2 have been used in Eqs. ()31II33|) . Note that in the limit — one 
finds uj^^acc — 1 for type of models. 

The predictions given by Eqs. ()31ll33p are shown in Fig. |3J One may be surprised that 
different models lead to the same results for most central collisions. This is because tUp'^^ 
goes to zero at Np"-' ~ A, as it follows from Fig. El (left). The predictions of the T-, M- 
and R-models differ because of their different response on the N'^^^ fluctuations. These 
fluctuations become small in the most central events. Therefore, the best way to study 
the mixing-transparency effects is the analysis of the centrality dependence of the particle 
number fluctuations in the projectile and target hemispheres. 

We now discuss an effect of a limited acceptance for the approaches with randomly 
fluctuating source rapidities. In this case the scaled variances in the projectile and target 
hemispheres are given by Eqs. ()10II13|1 . provided a width of a rapidity spectrum of particles 
emitted from a single source is narrow. In a general case, when a rapidity width of the 
source particles and an acceptance window are comparable in size, it is difficult to make 
analytical estimates. The problem can be solved in the limit of very narrow sources (a 
source width, A?/^,, is much smaller than the experimental acceptance interval. Ay). The 
hadrons created by "narrow" sources have correlated rapidities, but Eq. can be used 
for the scaled variances of the number of sources assuming that the source rapidities in the 
projectile hemisphere are not correlated. In this case our Eqs. fll0lll3j) yield: 

u^.^UT; rr) ^ 1, (34) 

1 1 



^-,acc{M; rr) ~ 1 + n. 



+ ■ ( ^ + 7 uj'P 



2.6 + 0.2 uj'P , (35) 



,2 4 

ujP_^acc{R; rr) ~ 1 + n_ [l - g* + g* ■ w^"''^] ~ 2.2 + 0.8 oo'^'' . (36) 



As before, we use ~ 0.4 and n_ ~ 2 in Eqs. ()34II36|) . The corresponding curves are 
plotted in Fig. El 

13 




FIG. 4: The dependence of the scaled variance of negatively charged particle multiplicity in Pb+Pb 
collisions at 158A GeV on the number of projectile participants, N^°'' , in the NA49 acceptance 



located in the projectile hemisphere. The Pb+Pb data 
comparison the result for p+p interactions at 158 GeV 
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111 are indicated by filled circles. For a 



131 is shown by the filled square. The 



displayed errors correspond to the sum of systematic and statistical uncertainties. The dashed line 
shows a dependence predicted within T-models (Eqs. pij) and ()34() ). the solid lines correspond to 
the mixing(fr) models (lower line, Eq. (jSU), and the mixing(rr) models (upper line, Eq. (j^SJ). The 
lower dashed-dotted line corresponds to the reflection(fr) models (|33j) ). whereas the reflection (rr) 
models are indicated by the upper dashed-dotted line. 

The experimental points for Pb+Pb collisions at 158 A GeV clearly exclude transparency 
and reflection approaches discussed here. The mixing model with random fluctuations of 
a source rapidity and a narrow width also strongly disagree with the data. A reasonable 
agreement is observed only for the mixing-hydrodynamical and WNM models. We remind 
that a large degree of mixing was found previously in the analysis of the pseudo-rapidity 



spectra of charged hadrons produced in d+Au interactions at ^/snn = 200 GeV 
WNM 0- 
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Note however that the WNM model, used here as a simple example to illustrate the MinF 
method, can not reproduce many observables connected to collective behavior of matter cre- 
ated in high energy nuclear collisions, like radial and anisotropic flows and the strangeness 
enhancement. On the other hand, these effects are at least qualitatively described by sta- 
tistical and hydrodynamical approaches. 

7. At the end several comments are appropriate. We considered three limiting behaviors 
of nuclear flows: transparency, mixing and reflection. In general, all intermediate cases are 
possible and they can be characterized by an additional parameter. Eq. introduces a 
mixing parameter a for the net baryon number with limiting cases Q;(T) = 0, a(M) = 1/2, 
and a{K) = 1 for T-, M- and R- models, respectively. Within WNM a parameter, a -(A?/) 
fl^ . defined for each particle species, i, and for each rapidity interval. Ay, was suggested. 
The limiting cases are again: aj{Ay;T) = 0, aj{Ay; M) = 0.5, and al{Ay; R) = 1. The 
values of the mixing parameter can be extracted by fitting the experimental data. 

The fluctuations of the participant number lead to the fluctuations of the center of mass 
rapidity {Ay ~ -l/2\og{Np 'yNP^"^). This alone may result in additional multiplicity fluc- 
tuations. We estimated that for the NA49 data discussed above the corresponding increase 
of the scaled variance is smaller than 5%. 

The MinF-method can be used independently of the degrees of freedom relevant at an 
early stage of collisions (e.g. hadrons at a low collision energy or quark and gluons at a 
high energy). This is because the concepts of the spectators and the participants as well 
as hadron multiplicity fluctuations are valid at all relativistic energies and for all collision 
scenarios. In the case of collisions of non-identical nuclei (different baryon numbers and/or 
electnc change to ba.yo„ .atios) one can t.ace flows^ot the conserved chatges by looMng at 
their inclusive final state distributions (see e.g. 20]). An interesting information can be 
extracted from collisions of two nuclei with different atomic numbers (see jisl]). 

In the case of identical nuclei only the MinF method can be used. It gives a unique 
possibility to investigate the flows of both the net baryon number and particle production 
sources. 

8. In summary, a method which allows to find out what happens with the initial flows in 
high energy nucleus-nucleus collisions was proposed. First, the projectile and target initial 
flows are marked in fluctuations (the MinF-method) in the number of colliding nucleons. 
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This can be achieved by a selection of colhsions with a fixed number of projectile participants 
but a fiuctuating number of target participants. This case is considered in details in the 
present study. Other selections are also possible. Secondly, the projectile and target related 
matter in the final state of collisions are distinguished by an analysis of fiuctuations of 
extensive quantities. We apply this method to the NA49 data on multiplicity fiuctuations 
of negatively charged hadrons produced in Pb+Pb collisions at 158 A GeV [11|. The results 
are consistent with the model which assumes a significant degree of mixing of the projectile 
and target fiows at the early stage of collisions followed by the hydrodynamical expansion 
and freeze-out. 
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